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ABSTRACT 

We have measured transit times for HD 189733b passing in front of its bright 
(V = 7.67) chromospherically active and spotted parent star. Nearly continuous 
broadband optical photometry of this system was obtained with the MOST (Mi- 
crovariability & Oscillations of STars) space telescope during 21 days in August 
2006, monitoring 10 consecutive transits. We have used these data to search 
for deviations from a constant orbital period which can indicate the presence of 
additional planets in the system that are as yet undetected by Doppler searches. 
There are no transit timing variations above the level of ±45 s, ruling out super- 
Earths (of masses 1 — 4M e ) in the 1:2 and 2:3 inner resonances and planets of 20 
Mq in the 2:1 outer resonance of the known planet. We also discuss complications 
in measuring transit times for a planet that transits an active star with large star 
spots, and how the transits can help constrain and test spot models. This has 
implications for the large number of such systems expected to be discovered by 
the C0R0T and Kepler missions. 

Subject headings: stars: individual (HD 189733) - planetary systems - methods: 
photometry - data analysis 



1 Bascd on data from the MOST satellite, a Canadian Space Agency mission, jointly operated by Dynacon 
Inc., the University of Toronto Institute for Aerospace Studies and the University of British Columbia, with 
the assistance of the University of Vienna. 
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Introduction 



While groundbased radial velocity (RV) and photometric transit surveys have unearthed 
more than 200 extrasolar planets in just over the last decade^], the ability to detect planets 
similar to the Earth, in size and mass, has so far remained out of reach of both of these 
methods. However, the potential to discover planets comparable to the mass of the Earth 
currently exists through measurements of transit timing variations (TTVs) in known tran - 
siting planetary systems (jAgol et al.ll2005l ; iHolman fc Murray! 120051 ; iHeyl fc Gladmanll2006l ). 
Additionally, nearly continous photometry from space currently offers the precision and time 
coverage to search for Earth-size transiting planets (e.g., Barge et al. 2005, Croll et al. 2007). 

A summary of the TTV method is as follows. The motion of a transiting planet whose 
orbit is perturbed by another planet in the system will not have a constant period. These 
chang es in period occur on the level of seconds to hours (jAgol et al.ll2005l ; IHolman fc Murray 
20051 ) depending on the mass and orbital parameters of the pertubing planet. By measuring 
the times for repeated transits of a known planet, one can at least constrain, if not unambigu- 
ously determine, the mass and semi-major axis of the perturbing body. Such TTV analyses 
can be performed on any exoplanet for which high quality photometry is available for a 
numbe r of transits. For the TrES-1 system an analysis of transit timings has been carried 
out by ISteffen fc Agoll (120051 ) with the result that any companion planet in an orbit nearby 
to the known close-in giant plane t must have a mass c omp arable to or less than th at of 
Earth. For the HD 209458 system, kgol fc Steffen] J2007h and iMiller-Ricci et al.l J2007h have 
placed limits on the existence of additional low-mass planets down to less than an Earth 
mass in certain resonant orbits, using data from the Hubble Space Telescope (HST) and 



the M icrovariability and Oscillations of STars (MOST) (IWalker et al.ll2003l ; iMatthews et al. 



2004 ) satellite, respectively. Here we apply a similar TTV analysis to the HD 189733 system 
to search for low-mass companion planets in orbits neighboring the known transiting planet. 

HD 18 9733 is currently t he brightest star (V = 7.67) known to harbor a transiting 
exoplanet (jBouchy et al.l 120051 ) . This fact along with its position on the sky and the short 
(2.2 d) period of its transiting planet makes it ideally suited for observations by the MOST 
satellite. MOST observed HD 189733 for 21 days in August 2006, monitoring 10 consecutive 
planetary transits. By determining the timings of these transits we are able to place limits 
on the presence of additional planets in this system down to a level of several Earth masses 
in certain orbits. 



Another characteristic of HD 189733 that makes it unique among the known transiting 



2 The Extrasolar Planets Encyclopedia: 



http: / /www.exoplanet.eu 
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systems is that the star has surface spots which modulate its optical light at a level of about 
3% during the stellar rotation cycle (Winn et al. 2006a, Matthews et al. in preparation). 
This level of variability is consistent with the fact that the star is known to be relatively 
active, with a chromospheric activity index S = 0.525 (IWright et al.ll2004l ). The presence 
of star spots on HD 189733 must be taken into account when fitting the light curves of 
the exoplanetary transits, since the shapes of individual transits can be affected (1) if they 
coincide with the planet's passage in front of a spot, or (2) if the underlying variability of 
the star is rapid enough to affect the fit to the transit. 

The CoRoT jBariinkood : barge et al.ll2005h and Kepler faorucki et al.lhooiliasri et al. 



20051 ) satellite missions (launched December 27, 2006 and scheduled for launch in February 



2009, respectively) are both expected to detect a large number of new transiting planets to 
add to an already growing list. One exciting prospect of these new systems is that each 
of them should be suited to the type of TTV analysis that we present in this paper, since 
high-quality observations of multiple consecutive transits will be available. Many of these 
newly discovered exoplanets may transit in front of chromospherically active stars like HD 
189733, since such stars are commonplace, and neither the CoRoT nor Kepler team has done 
a priori selection of quiet stars for their transit searches. In the following TTV analysis, we 
incorporate the intrinsic variability of HD 189733, in an attempt to quantify the effects that 
this may have on the measurement of transit times. This provides valuable experience for 
future observations of transits of active spotted stars by CoRoT, Kepler, and groundbased 
surveys. 



2. MOST Photometry 

HD 189733 was observed by MOST nearly continuously for 21 days from July 31 to 
August 21, 2006. The photometry was collected in MOST's Direct Imaging mode operating 
with a single CCD, where a defocused stellar image was recorded on a subraster of the MOST 
Science CCD. The exposure time for individual frames was 1.5 s (read time is negligible), 
and 14 consecutive images were "stacked" onboard the satellite to achieve high S/N with a 
time sampling interval of 21 s. The original plan for the duration of the HD 189733 run was 
14 days, and the first 14 days of observation have a duty cycle of 94%. When examination of 
the light curve showed the obvious spot modulation, the run was extended for 7 more days, 
by sharing each MOST satellite orbit (P = 101.4 min) with the next scheduled Primary 
Science Target. Therefore, the last week of data have a reduced duty cycle of 46%, and 
the observations were restricted to MOST orbital phases of highest scattered Earthshine 
(with a resultant increase in photometric scatter). Still, the time sampling and photometric 
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precision for the final week of observations remain excellent and are sufficient for analysis of 
the transit light curve. 



The data were reduced in the same way as described by iRowe et al.l (120061 ) for the 
transiting system HD 209458. The stellar fluxes were extracted by aperture photometry 
(aperture radius = 4 pixels; stellar image FWHM ~ 2.5 pixels), which gives better results 
than PSF (Point-Spread-Function) fitting. The raw instrument al light curve was decorrelated 
against the sky background as described by IRowe et al.l (120061 ) and also against the location 
of the PSF centroid on the CCD. 

The reduced light curve of HD 189733 is shown in the top panel of Figure [H It is 
immediately apparent that the star itself was variable during the MOST observations at a 
level of approximately 30 mmag. The shape and timescale of these variations are consistent 
with rotational modulation due to spots on the surface of HD 189733. 



3. Fitting the Transit Light Curve 

To fit and time the observed transits properly, the underlying variability of the host star 
must be subtracted from the light curve. We have adopted two independent approaches: (1) 
We smooth the out-of-transit (OOT) light curve using as a kernel a weighted average of 
the data covering one MOST satellite orbital period (about 101 min) centered on each data 
point. For the in-transit data, we linearly interpolate the smoothed data through the time 
of transit. Finally, we subtract the resulting smoothed light curve from the data. (2) We 
apply a filtered Discrete Fourier Transform to the data, in which we remove all power except 
at the orbital period of the exopl a net H P 189733b and its harmonics. We use the orbital 



period determined by lWinn et al.l (120071 ). 



The filtered light curves resulting from these two methods are shown in the middle and 
bottom panels of Figure [Q respectively. In the subsequent TTV analysis, we use both of the 
independently normalized light curves and find consistent results from each. 



3.1. Transit Model 



In computing the transit times, we compare th e MOST data against a model transit light 
curve, constructed using the formalism set forth in lMandel &: Agoll (120021 ) . To determine the 
system parameters describing this model — stellar mass and radius (M* and R*), planetary 
radius (R p ) , orbital inclination (i) , orbital period (P) , and stellar limb-darkening coefficients 
(ci — C4) — we procede as follows. 
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Since M* cann ot be determined from photometry alone, we adopt the value from 
Bouchy et al.l (120051 ) of M* = 0.82 ± O.O3M , which in turn was determined by fitting the 
star's observed spectral parameters against stellar evolution models. For i?*, R p , and i, we 
take advantage of the high-quality photometry provided by MOST to do an independent fit 
to these parameters. To determine the transit parameters for the HD 18973 3 system, we use 
a max imum likelihood analysis to fit the data with the analytic model of iMandel fc Agol 
(120021 ) and also add reflected light from the planetary companion (see Matthews et al. (in 
preparation) for details). For priors on the stellar and planetary mass and radi us as well as 



the or bital period and inclination of the planet, we adopt the values given in IWinn et al. 
(l2007h . We obtain stellar and planetary radii of 0.749 ± 0.009 R and 1.192 ± 0.019 R Jup 



(equatorial radius), respectively, and we find the orbital inclination angle to be 85.70°±0.11°. 

For the ab ove analysis, we fi xed the orbital period of HD 189733b at 2.2185733 ± 
0.0000019 days (IWinn et al.l 120071 1 . This value was determined by compiling all available 
transit times for fully observed transits over a period of more than a year, making it more 
precise than what could be obtained from a fit to the MOST light curve alone. In Section H] 
we present our own calculation of P by adding our transit times to the list of those previously 
available, and we find that the value we obtain is entirely consistent with the one we have 
chosen to use for our model transit light curve. 

Since the MOST telescope has a custom broadband filter covering the range 350 — 
700 nm, we cannot use standard tables of limb-darkening coefficients (such as those by 
Claretl (120001 )). which are only valid for standard filter systems. In addition, the quality 
of the MOST photometry is not high enough to observationally constrain limb-darkening 
parameters, especially given the wide bandpass of the instrument. We instead generate 
non-linear limb darkening parameters for HD 189733 from a synthetic s pectrum , calcu lated 
from model atmospheres using the ATLAS 9 and ATLAS 12 codes by Kuruczl (119951 ) and 
rewritten in FORTRAN 90 (J. Lester, private communication). We employ a model with 
T e ff = 5000 K, logg = 4.5, microturbulence of 1 km/s, and solar me tallic ity, consistent wit h 
the stellar parameters for HD 189733 given by lBouchy et al.l (120051 ) and iMelo et al.l (120061 ). 
The emergent spectrum is calculated at each of 17 values of fi = cos 9 from the center of 
the star out to the limb, where 9 is the angle relative to the normal. We then determine 
non-linear limb-darkening coefficients (c n ) as a function of wavelength according to the law: 



/(r) = J>(1 - M n/2 ) ■ 



n=l 



Finally, we multiply by the total throughput function for MOST, which includes the effects 
of the filter, optics, CCD, and its electronics. The resultant limb-darkening coefficients for 
the MOST bandpass, along with the other system parameters for the model light curve, are 
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listed in Table HJ In Figure [2], the transit model is plotted over the binned (normalized) 
data, which have been folded at the orbital period of HD 189733b. 



3.2. Transit Times 

Using the transit model described above, we find the best-fit center-of-transit times for 
9 of the 10 transits observed by MOST. Transit 7 is omitted due to the fact that MOST was 
switching observing modes at this time to the dual-star mode. For the 9 transits that have 
good time coverage, we compute the model light curve for times that correspond with each 
of the MOST data points. We next determine the center time for each transit at which the 
X 2 value for the fit to the data is minimized. 

The 1-cr error bars a re calculated us i ng a b ootstrapping Monte Carlo simulation similar 



to the one described in lAgol fc Steffenl (120071 ). according to the following sequence. For 



each transit we shift the residuals from the best-fitted transit model (and their associated 
errors) by a random number of points. We then add the new residuals onto the transit model 
and recalculate the center-of-transit time using the same procedure described above, thus 
maintaining the original point-to-point correlations. We perform this procedure 200 times for 
each transit, and the scatter in the obtained transit times sets our uncertainties. This allows 
us to account for the effects of correlated (red) noise in the MOST data set, which is due to 
the repeating pattern of stray light (satellite orbital modulation of scattered Earthshine) as 
well as to other systematics that remain in the light curve including any residual intrinsic 
stellar variations. 

The resulting transit times are listed in Table El We show these measurements in an 
O-C (Observed time - Computed time) diagram in Figure [3[ where the expected time of each 
transit (for constant orbital period) has been subtracted from its observed time to search for 
transit timing deviations. The expected times of transit are calculated from the period and 
ephemeris that we present in Section HI In Figure H] we once again present an O-C diagram 
of the transit times for HD 189733b but now including, as a reference, all available timing 
data from the literature in addition to our own. 

For the first 6 transits observed by MOST, the transit times we obtain are consistent 
with a constant orbital period at the 1-0 level. The final 3 transits show increased scatter 
and uncertainty due to the fact that the light curve does not fully sample those transits, and 
portions of either ingress or egress are not observed. Gaps in the light curve during times 
when the planet is passing over the limb of the star can have a large effect on the accuracy 
to which transit times can be determined. This is due to the fact that fits to partial transits 
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are highly sensitive to the values used for the limb-darkening parameters and the orbital 
inclination. If any of these deviate from t heir true values by even a small amount, the effect 
on the transit time can be substantial (see iMiller-Ricci et al.l (120071 )). We do not use the last 
three transits for our subsequent TTV analysis, as we expect the formal error bars in these 
cases to be underestimates. 



4. Transit Timing Variations in the HD 189733 System 

The transit timing data on HD 189733 from MOST show no variations on three scales: 
(a) no long-term change in P OT \ } in about 15 years at the 60 ms level; (b) no trend in 
transit timings during the first two weeks of the MOST run; and (c) no individual transit 
timing deviations at the 45 s level, with a string of 5 consecutive transits showing no timing 
deviations larger than 10 s. For illustration see Figures [3] and HI 



4.1. Long-term Variations in the Orbital Period of HD 189733b 



We first address point (a) from the list above. Long-term variations in the orbit al period 
could be caused by several eff ects, e.g., orbital decay of the planet HP 18 9733b (jSasselov 
20031 ) , or precession of its orbit (jMiralda-Escudal2002l ; iHeyl fc Gladmanll2006l ) . By combining 



archiv al Hipparcos photometry from 1991 a nd 19 93 with the discovery data of iBouchy et al. 



(120051 ). iHebrard fc Lecavelier Pes Etangsl (120061 ) have m easured the orbita l period of HP 



189733b over a 15-year baseline to be 2.218574±°;°™ d. IWinn et all (120071 ) independently 



-0.000010 

determined the orbital period for HP 189733b to be 2.2185733 ± 0.0000019 days. Their 
measurement was based on a compilation o f all available tran sit times for fully observed 
transits from their own work as well as from lBakos et al.l (120061 ). 



To further refine the value for the orbital period of HP 189733b and also to reduce the 
size of its associated error bar, we pres ent here a new calculation of the period. We follow 
the same method as IWinn et al.l (120071 ) and use all available timing data for full transits of 
HP 189733b, adding our own data points from the first six transits observed by MOST to 
the list. We fit the resulting 18 transit times with a linear function of the form 



T C (N) =T c (0)+PN 



(2) 



where T c (0) is the epoch and N is the transit number. We choose iV = as the fourth transit 
observed by MOST, since this is the transit time known to the highest precision. The orbital 
period and epoch obtained from the fit are 2.2185733 ± 0.0000014 d and JP 2453955.52478 
± 0.00010, respectively. Our new value for the period is entirely consistent with all previous 
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determinations. We have however, reduced the size of the 1-a error bar to 120 ms, now 
making this the most accurate value available for the orbital period of HD 189733b. 

The star HP 189733 is estimated t o have a mass of 0.82 M & with about 3% uncertainty 
( jBouchy et al.ll2005l ; iBakos et al.ll2006l ). With our stellar model for this mass and a planet 
mass of 1.15 Mj up , we can compute orbital decay rates . For the fastest possible orbital 
decay (high dissipation linear model (jZahnl 119771 . 119891 )). we get a rate of 0.6 msec per 
year. Currently available transit timing data for HD 189733b are not yet sensitive to long- 
term changes in P at this level. However, with the accumulation of future transit timing 
observations of this sytem, further constraints can be placed on this estimate. 



4.2. Short-term Variations in the Orbital Period of HD 189733b Due to 

Additional Close-in Planets in the System 



The existing radial velocity curve for HD 189733 (jWinn et al.ll2006l ) excludes additional 
periodic variations (other than the known giant planet) larger than 10 m/s in amplitude 
and up to a period of about 40 d. This puts a limit on the mass of a possible long-period 
perturbing planet at about M pcrt > 1O" 4 M ~ 0.lM Jup ~ 32M e . The TTV method is 
complimentary to the RV method in that it is generally more sensitive to smaller close-in 
perturbing planets. This is especially true in the case of the additional planet lying in a 
resonant orbit with the transiting planet, where the interactions between the two bodies are 
strongest. We use only the first 6 MOST transits here (see §3.2), so the run can constrain 
only short-period perturbers in nearby resonances to HD 189733b, and periods shorter than 
about 6 days. 

To determine what types of planets are ruled out by the MOST data, both in mass and 
orbits, we solve the classical N-body problem, 

d Xi Grrij(xi — Xj) 

— = -h j=1 .^ | ,3 , (3) 

where for 3 bodies, x describes the initial positions of the particles. A third body was inserted 
with an initially circular orbit, with periods ranging from 1 to 9 days in increments of 0.01 
days, with masses from 1 - 100 M e in 1 M e increments, and on a coplanar orbit to HD 
189733b. The solution was advanced at 1.0 second intervals for 100 orbits of HD 189733b 



2 x 10 7 s) using the LSODA routine from ODEPACK (IRadhakrishnan fc Hindmarsh 



19931 ). O-C values were then calculated by a linear interpolation to estimate the integration 



time when HD 189733b returns to the midpoint of crossing the disk of the star. 



To determine the limits that transit timing data can place on additional planets in the 
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HD 189733 system, we compute a Fourier Transform of the O-C series generated by the 
N-body code for each value of period and mass of the 3rd body. We then extract the largest 
amplitude as shown in Figure However, we note that the largest amplitude of TTV can 
only be recovered from transit timing data if the entire libration period of the two-planet 
system is fully sampled. Otherwise, there is a risk that the times of largest transit timing 
deviations could be missed or passed over by the observations. Figure [5] therefore shows 
what types of limits on companion planets could be achieved if data were available for the 
full libration period of each hypothetical two-planet system. Given that the MOST transit 
timing data only cover 6 orbits of HD 189733b, we use the following procedure to determine 
what additional planets are ruled out from the system for cases where the libration period 
is longer than the two- week duration of the MOST observing run. 

We compare the N-body results against the MOST transit timing data to determine 
the maximum mass that an additional planet in the system can possess, while still remain- 
ing consistent with the data. This is achieved by removing any linear trends from each of 
the O-C series generated by the N-body code, and then optimizing the agreement with the 
observed transit timing residuals by minimizing the chi-squared statistic. We then allow for 
the computed transit times to shift by an integer number of transits and determine the min- 
imum value for chi-squared once the overall best fit is achieved. For each configuration, the 
perturbing planet is ruled out if the minimum in chi-squared is higher than a set threshold, 
implying that the results from the N-body simulation are inconsistent with the MOST transit 
times. Figure M shows the results of this process. The maximum mass of perturbing planet 
that remains consistent with the transit timing data at a confidence level of better than 5% 
is plotted as a function of the orbital period of the perturbing body. This corresponds to 
a chi-squared of less than 11.07 for the best fit of the N-body simulation to the data, on 5 
degrees of freedom. 

According to Figure EJ the TTV limits exclude perturbers of greater than 20 M e for 
the exterior 2:1 mean motion resonance and greater than 8 M e in the 3:2 resonance. In 
the interior 1:2 and 2:3 resonances, additional planets with masses larger than 4 M e and 
1 M e are ruled out, respectively. A range of intermediate orbits near the transiting planet 
HD 189733b are unstable for the lowest-mass perturbers, resulting in ejections or collisions. 
Between 2.0 and 2.4-day periods we can rule out sub-Earth mass planets due to both TTV 
constraints and stability requirements. Additionally, we can place more stringent limits than 
the 32 M e value from radial velocity measurements in a number of non-resonant orbits with 
periods of up to 5 days. 

For a perturbing planet on an eccentric orbit, the TTV's induced in the orbit of HD 
189733b can be significantly larger than for the case of a circular orbit. Also, there is a 
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reduced range of stable configurations for planets on eccentric orbits, especially for orbits 
interior to that of HD 189733b. We have performed a limited number of additional N-body 
simulations for pertubing planets in both mutually inclined orbits relative to the transiting 
planet, and in initially eccentric orbits. In almost all of these cases, we can place even 
stronger limits on the presence of additional planets in the HD 189733 system. In this sense, 
the limits that we have placed on perturbing planets from Figure [6] are robust limits across 
the entire range of eccentricity parameter space, since we initially only considered planets 
with zero eccentricity. Additional planets residing in eccentric or inclined orbits would have 
even larger effects on the transit times of HD 189733 than what we reported above. 



5. Considerations Due to Star Spots 



The amplitude of the intrinsic stellar variations during the MOST observing run suggests 
that HD 189733 may have had fairly large star spots of high flux contrast during this time. 
It is therefore quite possible that the exoplanet may have transited one or more spots during 
the 21-day run. This type of event would alter the shape of the transit lightcurve in a 
predictible way, and we can look for te lltale signs o f such an occurrence. From deviations 
in the shape of the transit light curve, IPont et al.l (120071 ) observed HD 189733b occulting 
starspots in photometry from the HST. The approximate change in signal was 0.7%, and 
thus they inferred the presence of spots at least 80,000 km across (1.1 Rj up ), given an inferred 
temperature of the spots that was 1,000 K cooler than the star itself. 

Were the transiting planet to have crossed over a star spot during the MOST observing 
run, the depth of the transit at this time would be shallower than otherwise predicted due to 
the fact that the planet would be passing in front of a cooler, hence optically dimmer, region 
of the photosphere. Since the lessening in transit depth would occur only during the portion 
of the transit when the planet is actually in front of the star spot and not n ecessarily during 
the entire transit, the overall effect when fitting the iMandel fc Agoll (120021 ) transit model to 
the light curve would be twofold. First, the shape of the transit would deviate from the form 
predicted by the transit model, resulting in a poorer goodness-of-fit (x 2 ) value. Secondly, a 
shallower than normal overall depth would be needed to best describe a transit affected by 
a star spot. 

Taking into account the second effect, we have examined the MOST data for variations 
in depth between successive transits to determine if such a spot-crossing event took place 
during the MOST observing run. For each of the six fully observed transits, we now allow 
for the depth of the tr ansit, (i? p /i?*) 2 , to be a free parameter in our fit. While the standard 
Mandel fc Agoll (120021 ) template transit model cannot predict the correct transit shape in 
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the case where the planet crosses a star spot, this analysis still allows for an estimate as 
to whether there are significant changes in the average transit depth throughout the MOST 
observing run. To determine the size of the error bars of these measurements, we employ 
the same bootstrap method described above for the transit timing measurements. This once 
again allows us to account for correlated noise in the data, which can also affect the transit 
depth. The resulting values for transit depth are given in Table [3j The mean transit depth 
is 2.660% with a mean uncertainty of ±0.027%. All six measurements are consistent with a 
constant depth at the 2-a level. 

With the current precision in our transit depth measurements, we note that we cannot 
detect changes in depth smaller than about 3% at the 3-a level. Essentially, this means that 
the noise in the data could conceal the passage of HD 189733 over star spots smaller than 
approximately 0.2 Rj in radius, under the assumption that the star spots are 1,000 K cooler 
than the rest of the photosphere (see discussion below). It is entirely possible that this type 
of small spot is present on the surface of HD 189733. 



Croll et al.l (120071 ) have obtained the best fit of the unfiltered MOST light curve to a 



simple spot model with two large star spots, obtaining a rotation period of 11. 731^05 d 
and a stellar rotational inclination of 59lg degrees . The spot model per iod is very close to 
the rotation period of ll.953l0.09 d determined by lHenry fe Winn! (120081 ) from groundbased 
photometry spanning from October, 2005 to July, 2007. However, the two values still differ by 
several sigma; a fact that may point to differential rotation for the set of star spots that was 
being observed. The I Croll et al.l (120071 ) model for HD 189733 assumes circular spots, solid- 
body surface rotation, and the possibility of evolution in spot size during the observations. 
The star spots are assumed to be 1,000 K cooler than the rest of the photosphe re, which is 



consi stent with Sun spots and with HST observations of the HD 189733 system (IPont et al. 



20071 ). However, the star spot temperature cannot be determined uniquely from the MOST 
light curve, and there is a degeneracy between the sizes and temperatures of the spots. 

As a demonstration of the power of studying transits in systems with spotted stars, we 
combine the MOST transit observations with the star spot analysis of ICroll et al.l (120071 ) to 
predict whether HD 189733b should have transited in front of one of the modeled star spots 
during the observing run. The spot model predicts the position and size of the two large 
star spots on the disk of HD 189733 at any time during the MOST run. 

Using snapshots generated from the spot model discussed above, we can overlay the po- 
sition of the transiting planet on an image of the stellar disk to see if HD 189733b is expected 
to have passed over either of the two large star spots, as shown in Figure In this figure, 
it is possible to project the position of the planet on the star by the combination of three 
independent pieces of information. First, the sky projection angle between the orbital axis of 
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the transiting plan et and the stellar rotation axis of the star is known to be A = —1.4° ± 1.1° 
( IWinn et al.ll2006l ). This value was determined by measuring the Rossiter-McLaughlin ef- 
fect, which is the radial velocit y distortion seen a s a transit i ng planet occults the project ed 
rotation profile of the star fsee bhta et~aD koO$ \: iGimenezI (hood ): baudi fc Winnl (b()07t )). 
In addition, the stellar inclination angle (/*) has been obtaine d from combining the mea- 
surement of usm/* (also from the Rossiter-McLaugh lin effect) (IWinn et al.ll2006l ) with the 
rotation period of star obtained by ICroll et al.l (120071 ). Finally, the inclination angle of the 
planet's orbit is known from fits to the transit light curve. The only remaining ambiguity 
lies in the fact that it is not known whether the transit crosses along a chord below or above 
the stellar equator. In Figure H, we depict both cases. 

From Figure it can be seen that if HD 189733b were positioned above the stellar 



equato r as it transits, it would have crossed one of the large star spots in the ICroll et al. 



( 120071 ) model during 4 out of the 6 fully-observed transits. A complete occultation of a star 
spot comparably sized to the transiting planet would reduce the transit depth to 35% of 
its original value. Such a large effect would be readily noticeable by eye, and it has also 
been ruled out by the transit depth analysis presented above for the first 6 tran sits from 
the MOST light curve. We can therefore conclude that, if the ICroll et al.l (120071 ) model is 
correct, HD 189733b must transit below the stellar equator (as depicted in Figure [7j). In this 
case, the planet is not expected to cross over either of the large star spots, with the possible 
exception of Transit 3. In this single case, the combined uncertainties in the spot model and 
the orbital inclination axis of the planet may allow for HD 189733b to transit across the 
southern hemisphere spot. There is no evidence in the light curve of Transit 3 to suggest a 
full occultation of this star spot. However, a partial occultation cannot be ruled out at the 
precision of the MOST data, making t his transit the mos t likely one to be affected by the 
presence of star spots as determined by ICroll et al.l (120071 ). 



As an additional consideration, in the event that the transiting planet crosses in front 
of a spot during ingress or egress, the change in depth of the transit from its expected 
value would be minimal. However, since such an event would change the shape of ingress 
or egress, the effect on the measured transit time would be considerable as this is where 
the timing signal is most sensitive. This is not expected t o have occured in the HD 189733 
system during the MOST observing run if the ICroll et al.l (120071 ) spot model is correct (see 
Figure [7]). However, a transit across a star spot during ingress or egress could certainly take 
place during future observations. To determine the expected change to the transit time if 
such an event were to occur, we placed a simulated spot-crossing event into the data for one 
of the MOST transits. We found that a transit across a dark star spot with a diameter of 
0.5 Rj during ingress could alter the measured transit time by more than a minute. Such 
an event remains virtually undetected by eye in the MOST data, but could be detected with 
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spacebased photometry of sufficient time sampling and precision such as that of CoRoT, 
Kepler, or the HST. 



6. Discussion and Conclusions 

We have measured accurate transit times for 9 transits of HD 189733b observed by 
MOST during 21 days in August, 2006. For the 6 consecutive fully sampled transits we find 
no transit timing deviations above the 45-s level, allowing us to rule out certain configurations 
of additional perturbing planets in the system. We find no evidence for perturbing planets in 
the 1:2, 2:3, and 3:2 resonances down to the Super-Earth level (4, 1, and 8 M e respectively). 
In the exterior 2:1 resonance we rule out perturbing planets above the 20 M e level, and 
we place the most stringent limits available on a range of non-resonant orbits with periods 
ranging from 1 to 5 days. In addition, we find no evidence for long-term changes in the 
orbital period of HD 189733b between 1991 and 2006, indicating no significant observable 
orbital decay or precession. 

While HD 189733b is currently the only exoplanet known to transit a chromospherically 
active star, many more such systems may be discovered in the coming years, especially with 
the launches of the CoRoT and Kepler space missions. The occultation of a star spot changes 
the shape and depth of an observed transit, and passage of an exoplanet over a star spot 
during transit ingress or egress changes the measured center time of the transit. This is 
of particular concern for lower-precision photometry, where such deviations in the shape of 
the transit light curve may not be readily detectible, but the effects on the transit timing 
measurement can be significant. 

For the transit timing analysis of the MOST light curve for HD 189733, we have been 
able to account for the presence of star spots in several ways. (1) The nearly continuous 21- 
day coverage and high precision of the data make it possible to subtract the intrinsic stellar 
variability before performing the transit timing analysis. (2) We have test ed and constraine d 



a two-spot model for the longer-term variations in the MOST light curve (jCroll et al.l 120071 ) 



and combi ned it with spin-or bit alignment information from Rossiter-McLaughlin effect mea- 



surements (IWinn et al.ll2006l ). This has allowed us to project when and where the transiting 
planet would be expected to pass over a modeled star spot. (3) Our error analysis for the 
transit times takes into account correlated noise in the MOST light curve, due either to 
possible residual instrumental effects, or to passages of HD 189733b over small spots on the 
star that are below the detection threshold for this data set. 



Transiting planetary systems such as HD 189733 point to the intriguing possibility that 
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the transits themselves could be useful probes of the nature of spots and activity on the 
host star. High-quality photometry of future transits of this system could be used to 'map' 
the surface brightness and geometry of spots in the part of the stellar disk traversed by the 
transiting planet. 

We urge any observers reporting transit times for planets transiting active stars to 
carefully consider possible star spot effects in any observed deviant transit times before 
assuming the presence of a perturbing body in the system. The type of analysis presented in 
this paper can be implemented in principle for any system to rule out star spots as a cause 
of TTVs. This method necessitates a light curve taken with high photometric precision 
and complete time coverage over at least two stellar rotation periods, to detect and model 
the low-level variability of the star due to the motion of star spots across its surface. Such 
observations are best obtained from space or with a network of medium-aperture groundbased 



telesc opes such as the planned Las Cumbres Observatory Global Telescope (IBrown et al. 



20061 ). Such data will help to make headway in the search for low-mass planets in known 



transiting systems, even in the case where the stars are chromospherically active. 
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Table 1. Orbital and Physical Parameters for HD 189733 



Parameter 


Value 


M* [M ] 


0.82 ± 0.03 a 


R* [R©] 


0.749 ± 0.009 


R P z [Rjup] 


1.192 ± 0.019 


i[°] 


85.70 ± 0.11 


P [days] 


2.2185733 ± 0.0000019 b 


Cl 


0.825232 


c 2 


-0.977556 


c 3 


1.686737 




-0.657020 



From lBouchv et al.l (120051 ) 



b From lWinn et al.l (120071 ) 



Note. - C1-C4 are the non-linear 
limb-darkening coefficients calculated 
for the MOST bandpass. 
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Table 2. Best-Fit Transit Times 



Transit' ~n~ 


T„ (HTDl 

C 


a fHTDl 


-L I v~ 1. 1. L L ^. v V L X 


1 


2453948.86962 


±0.00043 


0.82 


2 


2453951.08753 


±0.00022 


0.83 


3 


2453953.30630 


±0.00030 


0.88 


4 


2453955.52467 


±0.00025 


0.61 


5 


2453957.74331 


±0.00027 


0.76 


6 


2453959.96201 


±0.00025 


0.65 


8 


2453964.39811 


±0.00044 


1.12 


9 


2453966.61827 


±0.00054 


1.00 


10 


2453968.83608 


±0.00037 


0.97 



Note. - Transit 7 has been omitted due to a large 
gap in the data at this time. 



Table 3. Transit Depth Measurements 



Transit # 


Depth (%) 


1 


2.661 ±0.048 


2 


2.605 ±0.015 


3 


2.631 ±0.019 


4 


2.670 ±0.028 


5 


2.701 ±0.023 


6 


2.694 ±0.021 
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Fig. 1. — Top: MOST light curve of HD 189733. Middle: light curve normalized by smooth- 
ing the OOT data and subtracting this smoothed curve. Bottom: light curve normalized by 
applying a filtered Fourier Transform to remove all power except at the orbital period of the 
planet and its harmonics. 
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Fig. 2.— Above: The phase diagram of the MOST photometry of HD 189733, and the 
transit model, folded at the orbital period of the planet. The data have been averaged in 
2-min bins. Below: Residuals from the model. 
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Fig. 3. — Deviation from predicted time of transit vs. transit number for transits of HD 
189733b observed by MOST in 2006. The first 6 transits were fully sampled, while transits 
8, 9, and 10 were only partially observed. 
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Fig. 4. — Deviation from predicted time of transit vs. transit number for all published well 
sampled tr ansits of HP 18973 3b. Data points from MOST are shown with filled symbols , 
those from I Winn et al.l (120071 ) with open red symbols, and those from iBakos et al.l (120061 ) 
with open blue symbols. Note that different methods have been used to compute the error 
bars for each individual data set. See the original papers for details. 
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Fig. 5. — N-body results for maximum transit timing deviation vs. orbital period of the 
perturbing planet. The color scale as defined on the right side of the plot indicates the mass 
of the perturbing planet (in M ffi ). The MOST data, covering 6 orbits of HD 189733b, show 
no transit timing variations above the level of 45 s. 
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Fig. 6. — Maximum mass allowed for a perturbing planet in the HD 189733 system, which 
still remains consistent with the MOST transit times. Planets occupying the region of 
parameter space above the curve are ruled out by the available transit timing data with 95% 
confidence. 




Fig. 7. — Snapshots from the HD 189733 star spot model by ICroll et al.l (120071 ). predicting 
the position of two large star spots at the times of the 6 complete transits observed by 
MOST. The projected tracks of the transiting planet are also shown, with two tracks for 
each transit since only the orbital inclination, not the direction (north vs. south) is known. 
A northern hemisphere transit is ruled out if the Croll et al. model is correct, since effects 
of the passage over the large star spot in this region would be readily visible in the transit 
light curve. 



